UUITP-14/09 



Spiky strings in the SL(2) Bethe Ansatz 

L. Freyhulf^Q, M. Kruczenski^^l and A. Tirziu^i 

" Department of Physcis and Astronomy, Uppsala University, 
P.O. Box 803, S-75108, Uppsala, Sweden 
'^Department of Physics, Purdue University, 
W. Lafayette, IN 47907-2036, USA. 

Abstract 

We study spiky strings in the context of the SL(2) Bethe ansatz equations. We find 
an asymmetric distribution of Bethe roots along one cut that determines the all loop 
anomalous dimension at leading and subleading orders in a large S expansion. At leading 
order in strong coupling (large A) we obtain that the energy of such states is given, in 

terms of the spin S and the number of spikes nhy E — S = ^In + 1^ „ 'n~ 

O (^^) . This result matches perfectly the same expansion obtained from the known spiky 
string classical solution. 

We then discuss a two cut spiky string Bethe root distribution at one-loop in the SL(2) 
Bethe ansatz. In this case we find a limit where n — > oo, keeping ^^f--, ^ fixed. 

This is the one loop version of a limit previously considered in the context of the string 
classical solutions in AdS^ x S^. In that case it was related to a string solution in the 
AdS pp-wave background. 
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1 Introduction 



Remarkable progress was achieved in the understanding of the AdS/CFT correspondence by 
studying the SL{2) sector of the theory. Much of this progress relates to twist two gauge 
theory operators of the type tr($D:^$). On the string side they are described by the folded 
string solution [1] as can also be verified by an alternative computation in terms of the cusp 
anomaly [2j. In the field theory a proposed all loop Bethe ansatz [3J can be used to compute 
their anomalous dimension at all loops in the planar approximation. Computations of quantum 
corrections to the folded string solution [1] at leading order in large S as well as its general- 
ization to AdS^ X O El [7] were recently used to successfully check such all loop asymptotic 
Bethe ansatz [HI [9l [TOl [11]. The all loop asymptotic Bethe ansatz computation was then ex- 
tended to the subleading orders in large S expansion [11], [121 [13] and precise matching with the 
corresponding string theory computations [H] was found. 

Beyond twist two, the SL{2) sector contains operators of higher twist which have also been 
investigated. For example operators of the type 



are described, on the string side, by the spiky string solutions [15]. There has been some recent 
progress in this area. In [TH] the solutions were shown to correspond to multi-soliton solutions of 
a generalized sinh-Gordon model. This allowed the construction of new, more general, solutions 
where the spikes move with respect to each other opening up the possibility to study a whole 
new class of operators. The relation with field theory was examined in detail in [17] where the 
elliptic curves associated with the classical solution where analyzed and a map was proposed 
to a similar structure emerging from the study of the field theory operators. 

However, with the current methods, full understanding of these operators from the field 
theory point of view requires the construction of the Bethe ansatz solution that describes 
them. It is the goal of this paper to fill this gap by providing a proposal to describe the spiky 
string solution in the all loop Bethe ansatz. 

It was shown in [HI [19] and extended in [17] that such operators can be described by a spin 
chain with a number of sites n being the same as the number of spikes. At leading order in 
large S the spiky string solution touches the boundary of AdS^. It was shown ^20j that at 
this order the spiky string solution can be mapped to the folded string solution by a conformal 
transformation, and thus the all loop energy at leading order in large S is given by 



where /(A) is a universal scaling function known as the cusp anomaly. In this paper we obtain 
this result from the field theory side by an all loop Bethe ansatz computation. 

In a semiclassical analysis in string theory the spiky string solution was extended to AdSs x 
in [21]. It turns out that, in the limit of large number of spikes n oo, there is an interesting 
scaling limit in which and - remain fixed. It was shown in |2T] that the solution in 

this limit is related to a solution in the AdS pp-wave background introduced in [2U]. In this 
paper, by using the found Bethe ansatz solution, it is shown that the same pp-wave like limit 
can be taken at weak coupling. 




(1.1) 
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To construct the Bethe ansatz solution, an important input is a set of integers Uk which 
may be interpreted as the bosonic mode numbers of the waves propagating in the spin chain. 
Different such sets of numbers give rise to different solutions of the Bethe ansatz equations. 
The folded string solution in AdS^ reduces [22] to the folded string solution in flat space in 
the limit ^ ^ 1. In fiat space we can quantize the string exactly in terms of right and 
left moving waves propagating on the string. It is easily seen that the folded string has the 
same number of right and left moving excitations, i.e. Ul = n^, with the same wave numbers 
kh = kR = 1. Extending this idea to AdS^, for all values of spin S and coupling A, suggests that 
in a Bethe ansatz computation the folded string solution is described by two sets of modes, 
one with = — 1 and the other with = 1. These numbers were used then as the input 
in the one loop SL(2) Bethe ansatz equations for two cuts (corresponding to Uk = ±1), and a 
symmetric Bethe root distribution was found to describe the folded string solution [23j. The 
same bosonic quantum numbers were used later [2H |3] to find the corresponding solution to 
the all loop Bethe ansatz. In the large S limit and finite J the two cuts in fact merge into one 
cut with a discontinuous wave-number distribution since we still have Uk = ±1. 

Following the same procedure, we look at the spiky string solution in fiat space, find the 
left and right wave numbers and use them as an input to find the corresponding Bethe root 
equation that we then proceed to solve. The spiky string solution was already discussed in 
[15j . It turns out that the bosonic quantum numbers that should be used in the Bethe ansatz 
equations are ni = —1 and = n — 1. Using this input we find an asymmetric distribution 
of Bethe roots spread over one cut along the real axis at one loop in weak coupling and leading 
order in the large 5* expansion. We then extend the computation to all loops and further obtain 
the all loop result also for the subleading order in large S. More precisely, we obtain that the 
energy of the spiky string at all loops can be written as 

E-S = n + fig)"^ fin 5 + 7i. + In f- sin -)) + "^B^ + n), / = (1-3) 

where -82(5') denotes the virtual scaling function of twist 2 operators obtained in [T2]. The 
leading large S term, i.e. E — S f{9)^ In S, matches the string theory result ( 11. 2p at the same 
order. 

If we further expand the all loop result (11. 3p for large coupling we obtain 

E-S = In^^ + ln -sm- -1 (1.4) 
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Remarkably, the first line precisely matches the result known from expanding the classical string 
solution |14j, that is, not only at order InS" but also S^. The second line is a prediction for 
the quantum corrections to the classical solution which will be interesting to check explicitly. 
Let us observe that, since the asymptotic Bethe ansatz reproduces the correct strong coupling 
result to order 5*°, wrapping effects play no role to this order in S. 

Having understood the solution with one cut we proceed to consider a two cut solution to the 
one-loop Bethe ansatz equations. Again we use the "spiky string" bosonic quantum numbers 
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for the integers rik, i.e. one cut has = —1 and the other = n — 1 where n is the number of 
spikes. The solution has four real parameters d < c < a < b which represent the position of the 
two cuts given by the segments [d, c] and [b, a] . The solution obtained allows us to compute the 
energy E, spin S, R-charge J and number of spikes n in terms of these parameters. Although 
there are four parameters it turns out that, for the solution to exist, a consistency condition has 
to be satisfied which implies that only three are independent. Equivalently this means that we 
can, implicitly, express the energy as a function of the other physical quantities E = E{S, J, n). 
As a check, in the limit c, 6 0, namely when the cuts merge, we recover the results of the 
previous one cut solution. 

More interestingly, in the limit d —>■ — oo a pp-wave type scaling is obtained, i.e. n —>■ oo, 
with and ^ remaining fixed. In this limit the equations simplify since one parameter 

is removed. To be able to compare with the string calculation of |2l] it would be necessary to 
extend this results to all loops and then consider the strong coupling limit. This should be an 
important test of the all-loop Bethe ansatz since these solutions contain a lot of structure and 
a highly non-trivial dependence of the energy with the other parameters. 

Another direction one can pursue is to consider the Bethe ansatz description of the more 
general solutions described in [TB|. In fact at this point it might be more interesting to find a 
direct mapping between the classical action of the string and the classical action of the spin 
chain (at strong coupling) in a similar way as in [25]. Mapping the actions would determine 
a mapping between the configurations of the string and those of the spin chain. In particular 
this would allow to understand how to see the spikes and their motion from the spin chain/field 
theory point of view. In the approach we follow here this cannot be seen since we map a 
semiclassical or coherent state on the string side to an eigenstate of the Hamiltonian (Bethe 
ansatz). A more detailed map would have to be constructed between coherent states on both 
sides. 

The paper is organized as follows. In section 2 we review the spiky string solution in flat 
space which leads us to the correct mode numbers to be used in the Bethe ansatz. In section 
3 we present in detail the construction of the one-loop Bethe root distribution for a one cut 
solution at leading order in large 5*. In section 4 we extend the computation to all loops using 
the all loop asymptotic Bethe ansatz. Section 5 deals with the extension of the all loop solution 
at leading order in large S to the subleading order 0{S^). In section 6 we analyze the spiky 
string solution with two cuts using the one loop Bethe ansatz equations. We then take two 
limits: first a limit when the two cuts collide recovering the one cut solution, and second, we 
consider the pp-wave type limit. In appendix A we present some details about the solution of 
the Bethe ansatz equations which we used throughout the paper. In Appendix B we present a 
review of the 0{S^) corrections to the ground state, which corresponds to the folded spinning 
string solution. 

2 Spiky strings in flat space 

As with the folded string, it is convenient to start by studying solutions in fiat space. In 
that case we can quantize the theory exactly and it turns out that the solutions are just a 
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superpositions of a left and a right moving waves: 

X = A cos((r2- 1) cr+) + y4(ra- 1) cos(cr_) (2.1) 
y = Asm{{n-1) a+) + A{n-1) sm{a^) (2.2) 
t = 2A{n-l)T = A{n-l){a+ + a^) (2.3) 

where = t + a, = t — a, n is the number of spikes and A is constant that determines 
the size of the string. Here, (r, a) parameterize the world-sheet of a string which is moving in 
a Minkowski space with metric: 

ds^ = -dt^ + dx^ + dy"^ (2.4) 

The solutions are periodic in a with period 27r and satisfy the equations of motion in confor- 
mal gauge, {d"^ — d1)X^ = 0, as well as the constraints (9+X)^ = {d^XY = 0. Quantum 
mechanically the state has riR = J^(n — 1)^ right moving excitations of wave number kn = 1 
and riL = A'^{n — 1) left moving excitations with wave number ki = n — 1 (satisfying the level 
matching condition rifikR = nikL)- 

All excitations carry one unit of angular momentum and therefore the total angular momen- 
tum and energy are given by 



S = nL + nR = A^n{n-l), E = ^2{nLkL + nnkR) = 2A{n - 1), E = 2^——S (2.5) 

which agrees with a classical computation. For n = 2 we recover the standard Regge trajectory 
E = \/2S and for n > 2 we get a Regge trajectory of modified slope. 

At fixed time, the shape of the string for different values of n takes the form depicted in 
figs. ([I]) and ([2]). It can be seen that the string has n spikes or cusps and, analyzing the time 
dependence, that it rotates rigidly in such a way that the end points of the spikes move at the 
speed of light. More details on these solutions can be found in [15j. Here, we will just need the 
fact that the spiky string is obtained by superposing left and right moving waves of different 
momentum. This suggests that when extending to all loops we should take a distributions of 
Bethe roots such that the wave number can take the values = — 1 or n^^ = n — 1 where n is 
the number of spikes. 



3 Spiky strings in the Bethe Ansatz: 1-cut solution 

We follow the same prescription used in [2l] to obtain the InS* dependence from the Bethe 
ansatz at large S. In the SL(2) sector the Bethe roots are on the real axis and, for large S", they 
accumulate on various segments or cuts. In the case of the folded string with J = 2 and S" 1, 
there is only one cut and the root distribution is symmetric^. In this case, namely for the lowest 
twist, there is only one state whereas for J > 2 there is more then one state. Nevertheless, for 
the lowest energy state at higher twist, the root distribution is again symmetric as discussed 

^In the case when S ^ J ^ 1 there are two cuts. This situation was analyzed in [5^ and will be studied 
later in this paper. 
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Figure 1: Rotating string with 3 spikes in flat space. We plot {x,y) in eq.( ]2.3p parametrically 
as a function of cr = — >^ 27r for n = 3, A = 1 and r = 0. 



in [21]. In this section we consider excited energy states which describe the n-spike solution in 
the large S limit for a fixed finite J. 

Following the discussion of the previous section, in the case of the spiky strings at leading 
order in large S, we assume to have one cut but with an asymmetric root distribution. The 
mode numbers of each excitation will be taken to be = — 1 or = n — 1 where n is the 
number of spikes. 

Let us start with the one loop SL(2) spin chain. The one loop Bethe ansatz equations 
corresponding to a XXX_i nearest neighbor spin chain can be written as 

u^y^^u^-u^^ *=1,2,...,S (3.1) 

Taking the logarithm of (13. ip . and further expanding in ^ as appropriate for large S gives [21] 

7 ^ 1 

— = 27rnfc - 2 V (3.2) 



where are the Bethe roots and J is the length of the corresponding spin chain. The total 
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Figure 2: Rotating string with 10 spikes in flat space. We plot {x,y) in eq.( ]2.3p parametrically 
as a function of cr = ^ 27r for n = 10, A = 1 and r = 0. 



momentum should vanish due to the cyclicity of the trace. This condition reads 

, i 



n 



The one loop energy is given by 



E-S~J 



(3.3) 



E;t-t (3.4) 



+ 

At this point one should solve eg. (13. 2 p for the real numbers Uk subject to the condition (13. 3p . 
However, when the number of roots S is very large, the roots accumulate in cuts on the real 
axis and one can conveniently approximate the problem by defining a density of roots. Thus we 
introduce the root distribution function po{u) = Xlfc"^!"^ — Uk)- Furthermore, for large S", the 
cut grows as 5* and we can make the approximation that 1. Therefore, at leading order in 
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Figure 3: Wave numbers n^('u) 
modes with -u < to have n,„ = 



for n=2,3,6. By analogy with fiat space, we take the 



—1 and those with u > to have n„ 



n — 1. 



the large 5" expansion (13.21) . we need to solve the equations 

Po{u) 



f 



du- 



u — u 



J dupo{u) = S 



du 



PoW 



0, 



u 



(3.5) 
(3.6) 
(3.7) 



which are, the Bethe equation, the normalization of the density, and the zero momentum 
condition. The function po{u) has support on a finite interval that needs to be determined as 
part of the solution. 

To solve these equations we use the method described in the appendix. We first define the 
cut, namely the interval where po is non- vanishing, to extend from d < x < a with d < and 
a > 0. The function n,,,(u) is defined as 



-1, d <u <0 
n — 1, < u < a 



(3. 



In figure 3 we plot such functions in terms of a rescaled u = ^. We then introduce a function 
F(w) = ^/w — a\Jw — d which has a cut on the real axis extending precisely from d to a. Fur- 
thermore, right above (below) the cut it is purely imaginary with positive (negative) imaginary 
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part. Using this function we define the resolvent as 



1 r 

G[W) = — — ^ 

^ ^ n F{u + tO+)u-w 

which, following the appendix A, determines po{u) to be 



du 



Po{u) = - 



Fiu) 



F u' 



du 

u' — u 



Besides, we need the behavior of G{w) at infinity which is readily seen to be 



GH -6-0 + ^1- + 
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Go 



riwiu) n — 2 in 

— —du = i 1 arcsin 

TT L F(u + zO+ 2 TT 



1 

TT 

a + d 



a + d 
a — d 



Go 



1 uriy^iu) 



-du 



2 " vr F{u + iQ- , 
As explained in the appendix we need to satisfy the consistency condition 

G{w) 



Res 



w — u 



oo 



-Go = 



for Pq{u) to solve eq.( 13.5p . On the other hand, equation ( ]3.6p gives 

S = po{u)du = f G{w)dw = ivrRes [G{w), oo] = —inGi 



and eq.f l3.7l) is 







Po[u 



du 

u 2 



1 / Giw] 



dw = ivrRes 



c 



w 



G{w) 



oo 



w 



(3.9) 



(3.10) 



(3.11) 

(3.12) 
(3.13) 



(3.14) 



(3.15) 



(3.16) 



which is satisfied if f l3.14p is. In the previous equations we used a contour C that encircles the 
cut in a clockwise direction. Using (13.121) we then find 



- = tan" ( —) 



d 



2nJ 



Using ( ]3.13p and ( 13.15^ gives 



S = nV—ad 

Before computing the energy let us compute the generating function G 



(3.17) 
(3.18) 



G{w) 



F(w) 



du 



= -z ± ^ In V^^+VE^^ (3.19) 

[u — w)^y {a — u){u — d) vr ^y —d{a — w) — ^y a{w — d) 
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where ± refer to w having positive or negative imaginary part, w takes values in the complex 
plane with a cut along the real axis from d to a. 

We now need to compute the energy which, within the Bethe ansatz is given by 



E-S-J 
2? 



d u^ + j 



G{w) 



m < Res 



G{w) i 



+ Res 



dw 



G{w) _i 



(3.20) 



^ { gC^) - G{-'-) )(3.21) 



Using the definition of G{w) and through an explicit computation we find 
E-S-J 



27rReG(^) = 2Im / I^'^^Mdu 
^2^ J \Fiu)\u-^ 



n In 



(^a(i -d) + 



;{a — d) 



(3.22) 

(3.23) 



For fixed n and large S we see from eqs. (13.171) and (13.181) that a~(i~S'^lsowe can 
approximate the last equation as 



E-S-J 
2? 



nln 



Sad 
a — d 



5>1 



(3.24) 



Therefore, through eqs. (l3.18p . (13.241) . we have solved the problem by writing E and S in terms 
of two parameters a and d related by the condition (I3.17p . Explicitly we can write 



d 

E-S-J 



S TT 

— tan — 

n 2n 

S TT 

cot — 

n 2n 

2n m I — sm — 

n n 



(3.25) 
(3.26) 
(3.27) 



Notice that up to now we have not needed the actual density po{u). However, in the following 
section it will be needed for extending the result to all loops and subleading terms in the large 
S expansion. Using eq. (l3.10p it can be written as 



Po{u) 



n 






In 


TT 






2n 


I 


n 




2n 




TV 



a/ a{u — d) + ^y—d{a — u) 



^,/a(v^^^^ — ^y—d{a — u) 



arctanh^/^, a> u>Q 



arctanh v^^" 



d\/a—u ' 



d<u<0 



(3.28) 
(3.29) 



Plots of the root distribution are shown in Figure 4. It is clear that S can be thought as 
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Figure 4: Root density po{u) = Po(f ) for = 2, n = 3 and n = 6. The larger the n the more 
skewed the distribution becomes. 



determining the scale of the variable u. For that reason it is convenient to define a rescaled 
variable u = ^. Defining new parameters 

a=— cot— , b = (3.30) 

n n n sm - 

n 

such that | = — (a-(-6), | = 6 — awe find another convenient way to write the density: 



n 



P-{u+ ~a)~a + _ (^ + 5)2^ 



Pq{u) = —\n. — ^^^^ (3.31) 



~ {u + a)a — \J}p — {u + a)'^\f\P' — 



Summarizing, the result from this section is that the energy of the spiky string at one-loop is 
given by 

= 2nln^ + 0(^°) (3.32) 

When n = 1 this indeed recovers the twist two result in |24j. The leading order expression is 
proportional to n similarly as at strong coupling. 
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4 Spiky strings in the all loop Bethe ansatz 

After understanding the root distribution at one-loop we now proceed to find a solution to the 
all- loop Bethe equations. We will follow the procedure used in [24j. The asymptotic all loop 
Bethe equations in the SL(2) sector read 

2Jarctan(2ufc) + zJln ^ = 2Tihk — 2 7 arctanf-u^ — (4.1) 

5/2 l_ 2 1 + - S/2 

where 9{u^ v) denotes the dressing phase and x^{u) are defined through m ± | = x^{u) + -^Sj^)- 
For the ground state the "fermionic" mode numbers are given by [21] 



n. = fc' + ^e(A;), A:' = ±1 ±| . . . , (4.2) 

where k' = k — e{k)/2. 

For the highest excited state there is no gap in the mode numbers [H] 



nk = k', k' = ±^,±^,...,±^^. (4.3) 



In the limit of large S we define the density of roots as 



du 

and the all loop equations can then be written as 



^ s'-m (4.4) 



m2 + 1/4 du 1 + c/V(x+(m))2 ' + 1 

f\ , , ,^ d r , l-^Va;+(M)a;-(u') ,A , 



At one-loop this reduces to 



2vrpo(n) - 2 f rf./ , ^"^.^ ^ . (4.6) 



Rescaling the roots, u = Su, we find in the large S limit 

2J , . /•^~" , . , . f 1 



^ 7r5(M) = 27rpo(n) " 2 y _ _ du'po{u) \^7r6{u - u') + ^- _ ) (4-7) 
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and the one-loop equation reduces to 



= 7iJ6{u) + -f[ " du'jp^^^. (4.8) 



— b—a 



[U — W 



For the highest excited state J is equal to the number of cusps, n, and we recognize the above 
equation as the derivative of the one-loop equation. 

To simplify the presentation of the all-loop computation we will omit the dressing phase and 
only restore its contribution in the end. We can now proceed with the all loop equation by 
splitting off the one-loop piece of the density, p{u) = po{u) + 2g'^cr{u), 



a(u' 



= 27ra(w) - 2 / 

J-oo [u - wy + 1 

+ 2?X '"^°^"^^l'^^^-Wx-(.')-(.)J 



+ / dv!a{v!)—\1%hi- ^' ^ ;. . (4.9) 



oo 



du \ 1 — g"^ / x^ {u')x (u) 



In the above we have extended the integral boundaries to ±oo in the terms containing the 
higher loop density, a{u). For the term containing the one-loop density we find, in the large S 
limit, 

1 f\ > / /^ ^ /^o-i I- 9Vx+iu)x-iu') 
' du'poiu)— (2i\n ' ^ ' ^ ' 



"^9^ Jd du \ 1 — g"^ /x~^{u')x (u) 

1 1 



^x^{u)x {u') X {u)x^{u') 
du \ x~^{u) x^{u) 



^^lm_Gi^ (4.10) 



Using fl3.19p we find 



and 



= 2'na{u) - 2 



G{±^) =±-\nS + ... (4.11) 
a{u') 



oo 



{u - u'y + 1 



^ ad f I 1 

+ nlnS^— — — - + 



du \x~^{u) X (u) 
r , , , „d f , I- g^/x+iu)x-{u')\ 

After Fourier transformatioij^ and a redefinition of the density, 

a{t) = -2ne*'^^a{t)\YiS, (4.13) 



^We use the convention f{t) — dw e ™^ f{w). Details on the Fourier transforms used here can be found 
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the result is 

ait) 

with the kernel defined by 



- (^K{2gt, 0) - ^ dt'K{2gt, 2gt')a{t')^ (4.14) 



K{t,t')=Ko{t,t') + Ki{t,t') (4.15) 
2 °° 

^o(t,0 = TT7$^(2n- l)J2„-i(t)J2„-i(0 (4.16) 



n=l 

2 °° 

^i^*^' '^') = 777 5Z 2nJ2„(t) J2n(t')- (4.17) 

n=l 

Including the dressing phase in the computation replaces the kernel above with 

K{t, f) = Ko{t, t') + K,{t, t') + Ka{t, t') (4.18) 

where 

K,{t,t') = 8/ / dt"K^{t,2gt")^^K,{2gt'\t'). (4.19) 
Further we find that the energy is given by 

= 8/na(0)ln5 + C(5°). (4.20) 

which gives 

E-S=-f{g)\nS + 0{S') (4.21) 
which matches the string theory result at leading order in large S |20j . 

5 Subleading corrections for the highest excited state 

The computation of the first subleading corrections for the highest excited state follow the 
same logic as the corresponding computation for the ground state but uses an asymmetric root 
density. The subleading corrections for the ground state were computed in [HI [12] using a 
method different than what will be convenient here. We will here use the same method as in 
[211 [3] but also include subleading corrections, for convenience we include the corresponding 
computation for the ground state in Appendix B. 

5.1 One loop 

The one-loop equation is given in (14.61) . where J should be set to n for the highest excited state. 
An approximate solution to this equation is given by fl3.28p which we expand for large S 

Po{u) = -ln^-^ln(u2) + -ln -sin- +C(1/^). (5.1) 



TT 27T IT \ n n 
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The above density captures only the leading order in the large spin expansion correctly. To get 
the subleading corrections of order we split off the above solution from the density, 

po{u) = poiu) +r{u), (5.2) 

and solve for r{u). We are now allowed to extend the limits of integration to ±00. 

/ 2 7r\ n 
= 27rr(M) + 2n In S* — n In(M^) + 2?T, In I — sin 



n n J + \/A 



-2 / dv[-\\iS In(n^) + - In - sm 



—00 



1x271 TT \n n J J {u — vY + 1 



00 



- 2 / dv "^^^ (5.3) 
The equation can be solved by Fourier transformation 

rit)=n(——^,--:]. (5.4) 



,l-e-l*l \t\J' 

Computing the energy including the first two orders in the large S expansion we find 
E-S-n _ r Mu) f-^^_ riu) 



2^2 + ^2 + 1/4 

2n In S* + 2n In ( — sin —\^n\ dt 



n nj \l-e-W \t\ 

2n ( ln^ + 7i; + ln ( ^sin^ j j . (5.5) 



5.2 All loops 

The all loop equation can be written as 



d x{i/2 + u) r pW) 

n— m — — = 2np[u) ' 



du x{i/2 — u) {u — u'y + 1 



+ 



Again we will omit the dressing phase, only to restore it in the end. Splitting off the one-loop 
density 

p(m) =Po(M) + 2(?2a(M), (5.7) 

where Pq{u) satisfies (14. 6p . we find 

n ( d ^ x(z/2 + n) ^\ ^ a{u') 

+ r du'd{u') — (21 In l^ZZ^iM^lM^ (5 8) 

7-00 du\ 1 - g'^/x+{u')x-{u)J 
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where 



1 — g'^/x~^{u)x {u') 



(5.9) 



1 — g"^ /x^{u')x {u) 

Note that the hmits of integration in the integrals containing the higher loop density can be 
extended to ±cxd to this level of approximation. To obtain the all loop energy, 

E-S-n 



29' 



p{u) 



x^[u] 



X (u) 



(5.10) 



we can follow a similar strategy as above. We split off the one-loop density and then extend 
the limits of integration in (15.101) to ±cxo. It is then clear that we only need the combination 
a{u)+a{—u) to obtain the energy. To write the integral equation for the symmetric combination 
(T(n) + cr(— n) we need the combination I{u) + I{—u) which we can split into two parts 



I{u) + I{-u) 



1 



+ 



du'po{u') 



d 



1 d 
2^du 
1 



2i In 



(u) x~{u) 
1 — g'^/x~^{u)x~{u') 
1 — g"^ / x^{u')x~{u) 



u''^ + l/Adu \x^{u) X (u) J 

Note that in the integral in (15. lip we have extended the integration limits to ±oo 
Fourier transformation of the equation (15. 8p gives 

^2vre-*/2 {M2gt) - 1) = vr (^(t) + ^(-t)) (1 - e"*) + ^ (/(t) + /(-t)' 



(5.11) 



^9' 

+A7rg^te 



-t/2 



dt' ( Ht') + a{-t') ) e-''/^K{2gt, 2gt' 



(5.12) 



where the kernel is given by (14.150 . For i(t) we find, to first two orders in the expansion, 
1 



I{t) + I{-t) 



Sim (^n 5 + 7ij + In sin ^ ^ te-*/'K{2gt, 0) 

{K{2gt,2gt')-K{2gt,0)). 



+ Sirte 



dt'- 



e - 1 

Using the result for the one- loop energy (15. 5p and introducing the density, a{t) 
we hence find for the all-loop equation 



(5.13) 



(5.14) 



a{t) 



t 



n 



1 

+ 2 



e* - 1 L 

oo 



K{2gt,0)- In 5 + 7s + In -sin 



TT 



n n 



n 



8gH 



{U2gt) - 1) 



dt'- 



n 



{K{2gt,2gt') - K{2gt,Q)) 



dt'K{2gt, 2gt')a{t') 



(5.15) 
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Taking the dressing function into account replaces the kernel by (14.181) . The string energy can 
now be written as 

E- S = n + lQg'^a{0). (5.16) 

We note that the integral equation is very similar to the equation in [12] and the solution can 
consequently be written as 

E-S = n + fig)"^ (^In ^ + 7^ + In sin ^ + ^^B^ig) + ... (5.17) 

where -62(5') denotes the virtual scaling function of twist 2 operators. Using the first orders in 
the expansions at strong coupling \12\ . 

f(g)=Ag-^ + Oil/g) (5.18) 
vr 

B2ig) = {-IE - \^g) fig) - - ln2) - f 1 - ^ + ^i^] + 0{l/g) (5.19) 



we find 



E — S = 2^(7 ( In — + In ( — sin — 

\ g \n n 

n 61n2 31n2, /4 . 7r\ 31n2, S\ , . 

+ - 1 + In -sm- In- +0(1/^) (5.20) 



2 \ TT TT \n n J 71 g 

We see that the leading strong coupling result is in agreement with the known string theory 
result IH]. For n = 2 we recover the result for the folded string [12J, this result is in full 
agreement with the one-loop computation from string theory [T^ |26] . For arbitrary n (15.201) 
is expected to be in agreement with the one-loop result computed from the sigma model [27j. 
Nevertheless, as we mentioned already, it would be interesting to compare the result in the 
second line in (15.201) with an explicit one-loop string computation of the 0{S^) term. 

6 Spiky strings in the Bethe Ansatz: 2-cuts solution 

We start again with equation (13.21) where we do not rescale u by S, i.e. we introduce the root 
distribution po{u) = J2k=i ^('^ ~ 



-f du \ = 7my,{u) - = 7cNy,{u) (6.1) 
J u — u 2u 



The condition that the total momentum should vanish gives 

u+ ' 



dupo(u)\n ^ = (6.2) 



The root density po is normalized as 

dupo{u) = S (6.3) 
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The 1-loop anomalous dimension is 

E-S-J 



du 



(6.4) 



Let us consider two cuts d < c < b < a. Extending the bosonic wave number distribution 
from the previous section we take 



Then the root distribution function is 



—1, d < u < c 
n — 1, h < u < a 



(6.5) 



_ 2n {c-u'){h-u') 1 



n-r^rib - a)(u' - d) . , , N :r^r - 

~ '^)n[77 -ITJ-, V' ^] + ~ a)K\r_ 



-{b-d){u' -a)' 



+ 



J 



-mu') - B,{u')) 



(6.6) 



for the left interval u' G [d, c], and 



Por{u) 



2n {u'-c){u'-h) 



+ 



TT y (m' - d){a - u') ^(^a-c){h-d) 



, (d — c)(u' — a) 
{d - a)Yi[) (f- -(, r] + {u' - d)K[r] 



(a — c){u' — d) 



for the right one u' G [6, a] . Here 



(6.7) 



(a — 6)(c — d) 



(a — c)(6 — (i) 

and we used the following definitions of the elliptic integrals 



K\m] 



da 



nfn, ml 



da 



1 



'0 vl — m sin a 
Also, we denote -Bi,i?2 the integrals 

B,{u') = ^\F{u')\A^, B^W) = ^\F{u')\A^ 

We compute the residue term in flA.12p and obtain 



l-msin^al-'^sin « 



(6.^ 



(6.9) 



(6.10) 



vrRes 



G{w) 



oo 



w — z 



z[nA^ + ^{Ai- A^)] +A2 + in- 1)^ 
Ztt 

a + b + c + d, . J . . 
2 ^ ' + 2^^^ ~ 



(6.11) 



where Ai integrals are defined in the Appendix A and can be written in terms of elliptic integrals. 
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As explained in Appendix A we need to take the above residue to zero in order for the root 
distribution (16.6116.71) to be a solution of (16.11) . This gives the following relationship 



nAi + ^iA^-A^) =0, 
Zn 



A2 + {n- 1)^3 = 



Then one can write J, n as functions of a,b,c,d 



J 



27inAi 



n = 1 



A, 



A5 - A4 A3 
The normalization condition can be computed using the function G as 



(6.12) 



(6.13) 



J po{u)du = iRes [G{w), 00] 



(6.14) 



nAi + -^(^4 - A5) +A2 + in- 1)^3 + -^(^2 - A3) + Aq + {n - 1)^7 



Taking into account equation (I6.12p the normalization condition reduces to 

^(A,-A3)+A, + {n-l)Aj = S 
The consistency of equation (16. ip gives the equation 



(6.15) 



J I du ^^ ^ = 27r / dupo{u)n{u) 



(6.16) 



The momentum condition (16. 2p expanded in large u gives 



du- 







u 



We use again the residues of the function G to compute the integral in (16.171) 



du 



Po[u) 



u 



ivrRes 



w 



G{w) 



00 



w 



+ ivrRes 



G{w) 







w 



(6.17) 



(6.18) 



We already computed the first term and set it to zero in (16. lip . What remains is to compute 
the second term in (I6.18p . We obtain 



du 



^ = VMi-A, -{n- 1)A, - ^As + ^A,) 



(6.19) 



Therefore the momentum condition (I6.17P gives the condition 

J 



\/ahcd[Ai + (n - 1)^5 + — (Ag - Ag)] = 

277 



(6.20) 
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All Ai integrals depend on the parameters a, b, c, d. Let us write the physical charges in terms 
of these parameters 

MA2-^Asf A2 E-S-J .pr^/^i ....X 

^ = ^^+ A,{A,-A,) -^'T,^ 2g^ = 2-I^^[^(2)] ^6.22) 

where we wrote the energy in terms of the function G whose explicit expression in terms of 
elliptic integrals is in flA.20p . 



We also have the condition (16.201) giving a relationship among parameters a,b,c,d 

^' - T,-"" + AM. -A,) = " 

We observe that using the above equations one can in principle solve for the unknown constants 
a,b,c,d and obtain the 1-loop anomalous dimension E = E{n,S,J). However, this is not 
possible to find explicitly because of the complicated elliptic integrals involved. Below we 
consider particular limits. 

6.1 Particular limit: 1 cut solution 

From the above general 2-cut solution for an arbitrary number of spikes let us recover the 1-cut 
solution again with arbitrary number of spikes. The 1-cut solution should be obtained in the 
limit 

6 = c = 0, d<0, a>0 (6.24) 

To obtain this limit we take b = e, c = —rje and take e — 0+. The parameter 77 > is to be 
determined. The integrals in this limit become 

Ine ^.on . TT a + d , vr a + d , , 

Ai = ^= + 0(e"), A2 = -- + arcsm -, A3 = - + arcsm (6.25) 

V —ad 2 a — d 2 a — d 

2cot"\/r7l , , 2tan~\/r7l , ^ , , 

A, = _V - + 0(ln6), A,= y - + 0(lne) (6.26) 

\/—ad^ e \/ —ad^ e 

I 1 ^ / 7\ + a + (i 

A^ = \J —ad (a + d) arcsm (6.27) 

4 2 Ci Qj 



. I : vr , (a + d) .a + d 

A^ = y —ad H — (a + d) arcsm 

4 2 d — a 



(6.28) 



. _ - - 1) cot 1 ^ 1 + (r] - 1) tan 1 ^ 1 

^8- 1== — — ^ + <^(inej, -= — — 

y—adriy/r] e y —adri^ri e 

Then the second equation in (I6.12p implies 



= — COS — (6.30) 

a — d n 
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which indeed matches the resuh from the section 3. The first equation in fl6.13p becomes 



J=— vmelne —-. _, 

tan y/r] + cot y/r] 

The normahzation equation (16.151) imphes 

^ n 7T 

b = [a + a) tan — 

2 n 







(6.31) 



(6.32) 



which again matches the results in section 3. Finally, the momentum condition (I6.20p deter- 
mines 7] to be 



n-2 1 

T] = 1 + 7T\I ^== + ... 



n vine 



(6.33) 



Equations (16.301 16.321) can be used to find the parameters a, d in terms of S*, n. The Bethe 
root density (16.61 16.71) in this limit reduces precisely to the 1-cut density expression (13.311) . 



6.2 Particular limit: "j3j9-wave" type scaling 

We want to consider the limit of large n. To achieve this we take the limit d 
keeping a, 6, c fixed. The integrals become 



2 1 



^Ja — c 



^2 = -TT + 0[ 



2 



A, 



-d 
2 



a — b 



a — c 



1 



..[cK[q] + {a-c)E[q]]^ + ... 

\/a — c V—" 



-^=[K[q] - n[-^, q]]-^ + o(^-^; 

a^a — c a — c \J —d [—dy^i^ 



2 ^.a — b, 1 



'-d 



A = ~d-^{a + b + c) + 0{-^: 



ay/ a — c a 

-- —^=[2{a -c){a + b + c)E[q] + (a(c - b) + c{b + 2c))ir[g]]-L= + 
3vct — c y —d 



= ^, [a(c - a)E\q] + b{a + c)ir[g] - {be + ab + ac)n[^— , g]]^i= + .. 
a^bcy/a — c a — c y —d 

Aq = ] [a(c - a)^[g] - acis:[g] + (be + ab + ac)n[- — q]]^= + ... 
a'^bcya — c a y —d 

The root distributions (16.6116.71) in the large —d limit reduce to 
, 2n 



(,.- - ^(n- - c) / a-^ ^ n[^.,]) + O(^) 

vrva — cV a — -u' \m' a — u' J y—d 



, .. 2n (u' — b)(u' — c) f a , ^,a — u' A ^/ 1 \ 

PoiiK) = , V -{-K[q] - n ,q] + Oi^=) 

nya — c V a — u' \u' a — c J \J —d 



-oo while 

(6.34) 

(6.35) 
(6.36) 
(6.37) 
(6.38) 
. (6.39) 
(6.40) 

(6.41) 
(6.42) 
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The number of spikes in this hmit is 



TT \/ a ~ c I — - 
n= — — ; r-zn— TV — « + ••• (6.43) 

2 cK\q\ + (a - c)E\q\ ^ ' 

In this hmit the physical parameters scale as 

n > 1, 5 ~ n^ J ~ E-S (6.44) 

which is the same as the pp-wave limit scaling found at strong coupling in [201 El]- More 
precisely we obtain 

^ n[!ti,9l + n|^.,l-A'|,l 

The condition (16.201) gives the relationship 

„,a — b , ab + be + ac , , , 

2n , q] = K[q] (6.46) 

a be 

The above three equations are to be solved for a, b, e as functions of S, J, n and then plug them 
in the energy. 

In the large d limit the function G{w) becomes 

X ^ za/(w - a){w - b){w ^c) H{w)^f^ 

w{w - a) ((a - e)E[q] + eK[q]){K[q] - n[^, q] - n[^, g]) 

(6.47) 

where 

H{w) = wK[qf + w{li[^^,q] + Ii[^,q])li[^^,q] (6.48) 

a a — e a — w 

+ K[q] [{a - w)U[^, q] + (a - w)U[^, + w)Il[^^, q] - aU[^^, q]] 

a a — e a — w a — e 

To find the energy we need to take w = ^ in G{w). The end of the right cut a is still large, i.e. 

3> a 3> 1, therefore we can expand in large a and finally obtain the energy at leading order 
in large n 

^ ^ = ^-^ (1 + 46^)^(1 + 4c^)^ sin[i(tan-i 1 + tan'^ h]{K[q] - n[^, q]) 

^9 va(a - c) 2 26 2c a 

(6.49) 

When b, e are also large, i.e \d\ 3> 6, c S> 1 the expression simplifies 

E — S — J nib + e) ,^^r. ^.a — b 

\/a — ey—abe a 

We therefore obtain a set of three equations fl6.45[ I6.46P which are to be solved for a, 6, c in 
terms of S, J,n. Furthermore, such solution for a,b,e should be replaced in the energy (16.491 
16301) . 
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7 Conclusions 



We have found the Bethe root distribution that describes the spiky strings in the context of the 
spin chain model that describes the field theory operators in the SL{2) sector. At one-loop we 
constructed sohitions where the roots condense on one cut and two cuts. The one cut solutions 
we were able to extend to all-loops. In particular expanding the result for large coupling we 
see that at order InS* and in the large S expansion, we recover exactly the result obtained 
from the classical solutions. The result contains a non-trivial dependence on the number of 
spikes and therefore is an important test of both, our proposal of describing the spiky strings 
by this particular Bethe ansatz solution and, more broadly, of the all-loop Bethe ansatz used 
in the calculation. In fact, the all-loop Bethe ansatz also provides a prediction for the one- loop 
quantum corrections to the spiky string. This is a doable calculation that should be interesting 
to perform in order to verify the prediction. In the case of the two cut solution, we did not 
extend the solution to all-loops. This should be interesting further work since the result at 
strong coupling is known from the string side. We have found that the equations simplify in 
a pp- wave-like limit where the number of spikes n grows to infinity keeping and ^ 

fixed. We showed that such a limit is well defined at one-loop STiggcsting that it actually can 
be taken also at higher loops. Finding the all-loop solution in this limit would be an interesting 
problem for further studies. 
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Appendix A: Solving the Bethe equations 

In this appendix we summarize the procedure that we used to solve the Bethe Ansatz equations. 
Although such methods are well known we adapt them here to our particular needs. However, 
it remains generic enough to be applied to other situations (e.g. at strong coupling.). 

To be concrete we have to find a function p{x) defined on the union of several segments 
Ci=i...n on the real axis (i.e. cuts) such that 

If Pj^^^^^^^ xeCiU...C„ (A.l) 

'"'Jciu...c„ X — X 

where j denotes principal part of the integral and N[x) is a given real function defined on the 
cuts. As a first step we have to find a complex function G{z), analytic except for cuts at Cj 
and a pole at infinity. Furthermore, we require that, on the cuts, the real part changes sign 
and the imaginary part is equal to N{x): 

G{x±i€) = ±p{x) + iN{x), xeCiU...Cn, e ^ 0+. (A.2) 
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The, yet undetermined, real part we denote as p{x). It is straightforward to find such a function 
in general for any N{x). While the procedure following below can also be directly extended for 
an arbitrary number of cuts, we take for simplicity two cuts. 

To be specific lets consider two cuts Ci = [d,c], C2 = [b,a] with d < c < b < a and the 
following function: 

F{w) = \/w — aV w — 6 V w — c\J w — d (A. 3) 

where the square roots are defined with a standard cut on the negative real axis. In figure [5] 
we show the phase of F{w) for w close to the real axis. Now we can write 



F{w) N{x) 
JciuC2 F{x + ie)x-w' 



dx 



(A.4) 



It is easy to verify that G{w) so defined has cuts in Ci and C2 and satisfies eq. (]A.2ll with 



p{x) 



7r 



F{x) 



F{x') 



N{x') 



■dx' 



X — X 



(A.5) 



Furthermore it behaves as G{w) ^ w{u} —>■ 00). Now consider an arbitrary point z (away from 
the cuts) and a small contour 7 encircling it as shown in figure [61 We have 



G{z) 

Deforming the contour we obtain 

G{z) = -Res 



2TTi F w — z 



dw 



Giw) 



w — z 



, 00 



1 

27ri 



G{w) 



■dw 



C1UC2 



w — z 



(A.6) 



(A.7) 



where the integral on the right hand side is over contours encircling the cuts clockwise. From 

'G{w) 



the property (1A.2I) we find 

G{z) = -Res 



00 



w — z 



+ 



JC1UC2 x' — z 



dx' 



(A.8) 



From this result and flA.2l) we find that, if z approaches one of the cuts from above, then: 



p{x)+iN{x) = G{x + ie) = -Res 

G{w 



—Res 



00 



w — X 



G{w) 
w — x' 

+ p{x) H 



00 



1 

+ — 



p{x') 



'C1UC2 



X — X — le 



-dx (A.9) 



C1UC2 



x' — X 



'Giw) 




, 00 




w — X 


^JCiuC2 



p(x') 



dx' , 



Therefore 

N{x) = iRes 

So we obtain that p{x) solves the problem ( lA.ll) under the condition 



dx', e ^ 0+ (A.IO) 



(A.ll) 



Res 



G{w) 



w — X 



00 







(A.12) 
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+1 ; — 4 1 — +1 

d +1 c b -1 a 



Figure 5: Phase of the function F{w) defined in the text, when w approaches the real axes. 
Along the two cuts, the imaginary part changes sign as depicted. 




Figure 6: Contours of integration surrounding z and the two cuts. 



The residue at infinity can be computed by expanding 



G{w) = Gow + Gi + -G2 + ... 

00 



which results in 



Res 



G{w) 



w — X 



, 00 



'Gq — xGi 



(A.13) 



(A.14) 



namely we need 6*0 = ^1 = 0. Since G{z) is already fixed we can only solve this problem for 
certain positions of the cuts. Namely, Gq = Gi = is an equation for a,b,c,d. In that case, 
formula flA.5|) gives the solution to the problem flA.ip . 

Furthermore, it is now straight forward to compute integrals of the type 



p{u) 
'C1UC2 ^ 
For example 



du 



Gjw) 

C1UC2 "y^" 



dw = vriRes 



G{w) 







+ TTzRes 



G{w) 



, 00 



(A.15) 



p{u)du 



C1UC2 



- ® G{w)dw = TciRes [G{w), 00] = —iTxG2 

2 JC^ VJCo. 



p{u) 
C1UC2 

p{u) 



du 



C1UC2 



u — z 



-du 



G{w) 
2 JciuC2 

G{w) 



C1UC2 



w — z 



dw = niG{0) — niGo 



dw = iTiG{z) 



(A.16) 
(A.17) 
(A.18) 



and so on. In particular the energy E and spin S can be computed by these formulas. Also, 
the last equation can be used as an alternative definition of G{z) |23]. In particular we can 
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recover flA.2l) by noticing that 



tnGiu'±ie) = [ ^ du = -f -f^du±inp{u') 

= -7iN{u') ±iTTp{u') (A.19) 



The above discussion of two cuts is vahd for any function N{u). However, in this paper N(u) 
is simply (16.11) N{u) = nw{u) — y~. For such N{u), the function G can be written in terms of 



eUiptic integrals as 



G{w) 



2i^J{w — a) {w — b){w — c) {w — d) 
7t{w — a){w — d)\J{a — c){b — d) 



(^2 



J 



2ttw' 



X ({a- d)n[^^^^^^^^-^, r] + {w- a)K[r]] - {n, ^ 



{b — d){w — a)' 

I / , -.^rid — c)(w — a) , , , 
X [{d-a)U[) -^,r] + {w-d)K[r] 



2txw' 



+ 



(a — c)(w — d) 
iJ \/ {w — a){w — b){w — c){w — d) 



(v44 - A 



(A.20) 



where ni, n2 are integers. In our case of interest rii = —1, n2 = n — 1. 



In the remaining of this appendix let us define the integrals used in this paper. First it is 
convenient to define the function 

\F{x)\ = y/\{x - a){x - b){x - c){x - d)\ (A.21) 

We can now compute 

dr 2 

A,= / -^ = ^= K[r] (A.22) 

Jt \F{x)\ ^^a-c){b-d) ^ ^ 

xdx 9 d — r 

A^= r^^ = ^ [(a-rf)n[ y ~''| ,r]+rfi^[r]] (A.25) 

A,= r -^^ = ^ ^ [(rf-a)n[ '^S^~'',| ,r]+air[r]] (A.26) 

/"'^ x^dx 1 
id ^{a-c){b-df ' 

d - c 

+ {a{a + c) + d{a-c))K[r]-{a-d){a + b + c + d)U[ ,r]] (A.27) 
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+ icia + c) + b{c-a))K[r]-{c-b){a + b + c + d)U[^^^,r]] (A.28) 



a — c 



As = / ^, = =\ad(a -c)(b- d)E[r] (A.29) 

J, x^\F{x)\ aWc^{a-c){b-dy ^ ^ ^ J K ) 

+ bd{—a^ + cd + a{c + (i))i^[r] + (a — d){bcd + acrf + a6c + a6(i)n[— | r]] 

+ 6(i(-c^ + ab + c{a + b))K[r] + (c - b){dab + ca6 + acrf + cbd)U[^ — ^,r]] 

Appendix B: Sub leading corrections for the ground state 

In the main text we compute the subleading corrections (in the large S expansion) to the 
density of roots. It is convenient to remind ourselves how this is done in the case of the ground 
state which is described by a symmetric root distribution. We do this at one-loop and and then 
make us of that solution to construct the all-loop solution. 

One loop 

The one-loop equation for the ground state reads 

27rpo(u) + 2n{J - 2)6{u) - - 2 dv P'^""] = 0. (B.l) 

+ 1/4 J_si2 [u -vY + 1 



An approximate solution to the equation is given by the following density [28] [21] which we 
expand for large S 



Mu) = - In l + ^l ^flll = ^ In 5 - i Hu^) + 0{l/S% (B.2) 

This density captures the leading order in the large spin expansion correctly. To get the 
subleading corrections of order we split the density, 

po{u) = pq{u) + r{u), (B.3) 

and solve for r{u). We are now allowed to extend the limits of integration to ±00. 

27rr (m) + 4 In S - 2 \n{u'^) + 2n{J- 2)5{u] 



00 



00 



m2 + 1/4 

-2 1 dv{-\nS-- \n(v^)] \. 2 dv- "-^^ . (B.4) 
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The equation can now be solved by Fourier transformation 

e-l*l/2 



f{t) = 

Computing the energy we find 
E-S-J 



2g' 



■S/2 + 1/4 



1 



+ 



(B.5) 



du 



r[u) 



41n^ + 41n2 + 2 



■oo 
oo 



dt 



m2 + 1/4 



1*1/2 



l-e-l*l 



\t\ 



-iJ-2) 



dt 



Ht|/2 



4(ln5 + 7£; - (J - 2)ln2) 



1 + e-l*l/2 

(B.6) 



All loops 

The all loop equation reads 



s/2 

■S/2 {U - m')^ + 1 



+ 



Ls/2 ^du\ 1- 



1 — g'^/x^{u)x {u') 
g'^ / {u')x~ {u) 



-2e{u,u] 



(B.7) 



In the following the dressing phase will be omitted for simplicity. It is easily restored in the 
end. Splitting the density 

p(m) =Po(M) + 2/a(M), (B.8) 

where po{u) satisfies fIB.ll) we find 

J . ^ . ^ . . _ . _ . ^^^^^ 



d , x(i/2 + u) 
— In 



2^2 ydu x{i/2-u) v? + l/A 
du \ 



2txo 



{u - v!f + 1 



, / 1 - g'^/x'^{u)x (n') 

+ / duaiu)— 2dn ' ^ ' ^ ' 



where 



1 f^'^ d f 1- 



1 — g"^ /x^{u')x (u) 

1 — g"^ / x^ {u)x^ {u') 
g'^/x~^{u')x~{u) 



(B.9) 



(B.IO) 



Note that the limits of integration in the integrals containing the higher loop density can 
extended to ±oo as we are considering the large S limit. Further we can split the integral, 
J(n), into two parts 



I{u) 



TC 



1 



+ 



du Pq{u') 



d 



1 d 
2^d^ 
1 



2i In 



(u) x^{u 
1 — g^/x~^{u)x~{u') 
1 — g'^/x~^{u')x~{u) 



u''^ + l/Adu \x^(u) X (u) 



(B.ll) 
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where G{u) is the one-loop resolvent, defined similarly to fl3.9p . Here we have used that the 
root distribution is symmetric around the origin. Note that in the integral in flB.ll|) we have 
extended the integration limits to ±00. 

Fourier transformation of the equation fIB.Qp gives 

±2ne-'/^M2gt) - 1) = 27r^(t)(l - e"*) + J(t) 

POO 

+ 8ngHe~'/^ / dt'a{t')e-''^^K{2gt,2gt'), (B.12) 
Jo 

with the kernel given by (14.151) . For /(t) we find, to first two orders in the expansion, 

i{t) = 16 (In 5 + 7b - {J -2)\n2)TTte~'/^K{2gt,0) (B.13) 
+ 4vrte-*/2 dt' (j-^ - (^(2^^, 2gt') - K{2gt, 0)) . 

Redefining the density as 

a{t) = -\e~'l'h{t) (B.14) 
we hence find for the all-loop equation 

a(t) = -±^\K{2gt,0){\nS + 'jE-{J-2)ln2)--^{Jo{2gt)-l) 
e* — 1 L SgH 

1 r°° / 2 J — 2 \ 

+ 2/ dt'[-^^~^^^)iKi2gt,2gt')-Ki2gt,0)) 

POO 

- Ag^ / dt'K{2gt,2gt')a{t') . (B.15) 







The energy is then given by 

E - S = J + lQg^a{0). (B.16) 

The effect of the dressing phase can now be restored by replacing the kernel above with (14.181) 
and we have hence obtained the same integral equation as in [12] . 

Note that this derivation explicitly made use of the fact that the density is symmetric. In 
the main text we do a similar calculation for the asymmetric distribution we are considering in 
this paper. 
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